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Abstract. 

We summarize the analytic model and numerical simulations of stochastically forced planets 
in a turbulent disk presented in a recent paper by Rein & Papaloizou [4]. We identify two modes 
of libration in systems with planets in mean motion resonance which react differently to random 
forces. The slow mode, which mostly corresponds to motion of the angle between the apsidal lines 
of the two planets, is converted to circulation more readily than the fast mode which is associated 
with oscillations of the semi-major axes. 

We therefore conclude that stochastic forcing due to disk turbulence may have played an impor- 
tant role in shaping the configuration of observed systems in mean motion resonance. For example, 
it naturally provides a mechanism for accounting for the HD12831 1 system for which the fast mode 
librates and the slow mode does not. 
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1. INTRODUCTION 

Of the recently discovered 336 extrasolar planets, at least 75 are in multiple planet 
systems [6]. More than 10% of these are in or very close to a resonant configuration 
where two planets show a mean motion commensurability, with at least five systems in 
or near a 2: 1 resonance [7]. 

Resonant configurations can be established by dissipative forces acting on the planets 
which lead to convergent migration [2]. An anomalous effective kinematic viscosity V, 
but with considerable uncertainty, has been inferred from observations of accretion rates 
onto the central protostars. The magneto rotational instability (MRI) is thought to be 
responsible for this anomalous value of v but the level of MRI turbulence and associated 
density fluctuations are uncertain. The density fluctuations give rise to random forces on 
close-by protoplanets. 

The influence of stochastic forces on migrating planets was explored first by [3]. 
They considered MRI simulations directly and therefore the simulation ran only for a 
relatively small number of orbits. Recently [4, and see references therein] presented 
a parameterized analytic model which can be used to predict the lifetime of resonant 
systems as a function of the stochastic diffusion coefficient. We shortly summarize this 
model in section 2 and focus on the stability of the GJ876 system in section 3 and the 
formation of the HD 12831 1 system in section 4. We refer the interested reader to [4] for 
a more detailed discussion. 



2. ANALYTIC MODEL 



Our goal is to estimate the growth rate of orbital parameters and the lifetime of resonant 
systems analytically without making any restrictive assumptions about the orbits or 
the nature of the stochastic forces. We therefore add an additional term to the full 
Hamiltonian of the system that is linear in a stochastic force F. The new equations 
of motion can be calculated by elementary means and expressed in terms of the orbital 
parameters such as semi-major axis a, eccentricity e, periastron GJ, mean motion n, mean 
longitude A. Thus, the additional terms due to the stochastic forces are linear in F. 

We assume that each component of the force F on each planet, in cylindrical co- 
ordinates, satisfies the relation Fi(t)Fi(t') = (F 2 ) g(\t — where the autocorrelation 
function g{x) is such that / °° g(x)dx = t c , with x c being the correlation time and (i 7 - 2 ) 
the mean square value of the i component. The stochastic forces make quantities they 
act on undergo a random walk. The square of the change of such a quantity A occuring 
after a time interval t is given by (AA) 2 = D t t, where = 2(F 2 ) t c is the diffusion 
coefficient. 

Considering stochastic forcing on a single isolated planet when Dj = D is independent 
of i, we obtain precise statistical estimates for the growth of the orbital parameters 
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When two planets are involved there are two modes of oscillation (libration) which 
we separate and describe as fast and slow modes with frequencies 0)/y and ft)/ v . To 
summarize our results, the ultimate lifetime of a resonant configuration of planets is 
determined by the time taken for the fast angle to achieve circulation. We determine this 
lifetime using equations similar to those described above but for the libration amplitude, 
obtaining 



Here the first quantity in brackets represents the ratio of the square of the central force 
to the mean square stochastic force acting on the outer planet. The other quantities in 
brackets are expected to be of order unity, while the last factor q/qcj is the ratio of the 
planet mass to the stellar mass compared to the same quantity in the GJ876 system. P\ is 
the period, with the subscript 1 denoting the outer planet. We showed that this is in very 
good agreement with numerical results [4]. 



3. STABILITY OF GJ876 

We have performed numerical simulations of two planet systems that allow for the 
incorporation of additional stochastic forces with the properties described above. In 
order to mimic the effects of turbulence it is necessary to calibrate these forces with 
reference to MHD simulations. The only free parameters are the mean square value 
of the force components per unit mass (F^) and the auto correlation time x c . The 




FIGURE 1. Time evolution of the resonant angles <j)\ = 2X\ — X2 — ®2 (fast mode), (fa = 2X\ — X2 — CJi 
and £ = Q2 — di (slow mode), the period ratio P1/P2 and the eccentricity, e\ , in the GJ876 system. The 
vertical lines indicate when the angles enter circulation for a prolonged period. The realization illustrated 
in the right panel scatters shortly after the fast libration angle, §\ goes into circularization. 



numerical noise generator that we implemented, uses a discrete first order Markov 
process [1] to generate a correlated and continuous force. We conclude from numerical 
MHD simulations that the natural scale for the diffusion coefficient is specified through 
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Of course we emphasize that the value of this quantity is very uncertain. 

In order to illustrate the evolutionary sequence we plot results for two realizations 
of the GJ876 system in figure 1. For these runs we adopted a diffusion coefficient of 
D = 0.42cm 2 /s 3 . In this context we note that reducing D increases the evolutionary 
time which has been found, both analyticaly and numerically to be l/D. The times 
at which the transition from libration to circulation occurs are different for the slow and 
fast angles and indicated by vertical lines in figure 1 . The tendency for the occurrence of 
very small values of e\ can be associated with transitions to circulation of the slow mode, 
£ being the angle between pericenters. We have verified this by considering the results 
from simulations which started with a larger value of e\. As expected, the attainment of 
circulation of the slow angle takes longer in this case. Also as expected, the time when 
01 attains circulation is not affected by the change in e\. 



4. FORMATION OF HD128311 

The ideas presented above can help us to understand the orbital configuration of the 
planetary system HD 128311. This system is in a 2:1 mean motion resonance with the 
angle §\ librating and the angle £ circulating (no apsidal corotation) [8]. 

It can be shown [2] that the planets should exhibit apsidal corotation if the commen- 
surability was formed by convergent migration only. Accordingly, it has been suggested 
[5] that an additional perturbing event (a close encounter with an additional planet) could 
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FIGURE 2. The plots show two different formation scenarios of HD12831 1 including turbulence and 
migration. We plot the observed system on the right hand side of each plot. The plots show the resonant 
angles (j>i, fa and £, the eccentricities e\ and ei as well as the period ratio P1/P2. 



produce orbital parameters similar to the observed ones. This perturbation is needed to 
alter the behaviour of £, so that it undergoes circulation rather than libration. 

We showed above that stochastic forcing that results from turbulence driven by the 
MRI readily produces systems with commensurabilities without apsidal corotation. This 
suggests that a scenario which forms the resonance through disk induced migration 
might readily produce commensurable systems without apsidal corotation if stochas- 
tic forcing is included. We present such formation scenarios that do not invoke artificial 
perturbation events in figure 2. We indeed find that model systems with orbital parame- 
ters resembling the observed ones are formed without difficulty and we are able to obtain 
similar final states for a wide range of migration parameters. 

With radial velocity measurements continuously improving and dedicated telescopes 
such as Kepler many exoplanets will be discovered within the next few years. That will 
allow us to assess the role of stochastic forcing on protoplanets in a statistical manner. 
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